PLB-GG-8192 



March. 24, 1996 



Strong supersymmetric quantum effects on top 
quark production at the Fermilab Tevatron 

Chong Sheng Li a ' b ' e , Hong Yi Zhou a ' c , Yun Lun Zhu a ' 6,/ and Jin Min Yang a,d 



a CCAST(World Laboratory) 
P.O.Box 8730, Beijing 100080, P. R. China 

b Department of Physics , Peking University, 

Beijing 100871, PR.China 
c Department of Physics, Tsinghua University, 
Beijing 100084, PR.China 
d Physics Department, Henan Normal University, 
Xin Xiang, Henan 453002, P.R.China 
e E-mail: csli@pmphpsvr.pku.edu.cn 
* E-mail: zhuyl@sun.ihep.ac.cn 



ABSTRACT 

The supersymmetric QCD corrections to top quark pair production by qq anni- 
hilation in pp collisions are calculated in the minimal supersymmetric model. We 
consider effects of the mixing of the scalar top quarks on the corrections to the total 
ti production cross section at the Fermilab Tevatron. We found that such correction 
is less sensitive to squark mass and gluino mass than in no-mixing case, and in both 
cases the corrections can exceed 10% even if we consider the recent CDF limit on 
squark and gluino masses. 
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The top quark has been discovered by CDF(DO) collaboration at Tevatron[l]. The 
mass and production cross section are found to be, respectively, 176±8(stat) ± 10(syst) 
( y 199t 1 2 %stat)±22(syst))GeV and 6.8^;|(6.4±2.2)p6 by CDF(DO) collaboration. At the 
Tevatron,the dominant production mechanism for a heavy top quark is the QCD process 
qq — > tt[2] . Recently, there has been a lot of interest in the one-loop radiative corrections 
of the top quark production cross section at the Tevatron, which arise from the new 
physics beyond the SM such as the two-Higgs-doublet model(2HDM) and the minimal 
super symmetric model(MSSM)[4,5]. In Ref.[5], the supersymmetric QCD corrections 
to the top quark production in pp were calculated in the simplest case of unmixed 
squarks and degenerate masses. But we have to know the impact from the mixing 
squarks because of interest to phenomenology, especially, sensitive dependence of the 
supersymmetric QCD on the squark masses and gluino mass [5]. In that reference also 
only two cases of the gluino mass m g = 150 GeV and rrig = 200 GeV, respectively, was 
considered. The purpose of the present letter is to evaluate the supersymmetric QCD 
corrections to the top quark production at the Tevatron in the general case of the mixing 
squark masses and compare our results with that in the case of unmixed squark masses 
given in Ref . [5] . We also discuss further the dependence of such corrections on the gluino 
mass. 

In the MSSM the strong supersymmetric interaction Lagrangian relevant to our 
calculation is given, in the present of squark mixing, by [6] 

C m = -iV2g s T a (qP R q L - qP L q R )~g a + H.C., (1) 

where g s is the strong coupling constant, T a are SU (3) c generators, Pl,r = |(1 ±75)) 9a 
are the Majorana gluino fields, and (jL tR are the current eigenstate squarks, which are 
related to the corresponding mass eigenstates qi ;2 by 

<?i = Ql cos 6 + q R sm6, q 2 = —q L sin 6 + q R cos6 (2) 
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The mixing angle 6 as well as the masses 2 of the physical squarks are determined 
by the following mass matrices [7] 

M 2 _ / m h + m t + 0.35D Z -m t (A t + cot (3) 1 

*~\ -m t (A t + //cot/3) m? fl + m\ + 0.16D Z ) ' ^ 

where = M§ cos 2/3 with tan (3 being the ratio of Higgs vcauum expectation values, 
m t L ,tR are so ^ breaking masses, A 4 is soft breaking parameter describing the strength 
of trilinear scalar interactions, and ji is the supersymmetric Higgs mass parameter. The 
relevant Feynman diagrams contributive to one-loop supersymmetric QCD corrections 
to qq — > ti amplitiude are shown in Fig.l of Ref.[5]. In our calculation we will follow the 
same notation and adopt the same regularization and renormalization scheme as was 
used in Ref.[5]. Notice that the off-diagonal elements of the squark mass matrices are 
proportional to the quark mass. In the case of supersymmetric partners of light quarks 
mixing between the current eigenstates can therefore be neglected. So we will take into 
account only the mixing between the squarks for the supersymmetric QCD corrections 
to the gti coupling and neglect such mixing for the gqq couplings because these couplings 
mainly involved light quarks. The renormalized amlitude for qq — > ti can be written as 

M ren = M + 5M 1 vertex + 5M 2 vertex + 5M hox (4) 

where M is the amplitude at tree-level and SM represent the supersymmetric QCD 
corrections arising from the effective gqq(gti) vertex and box diagrams, which are given 
by 

M = ^ P2 )(-tg s T a YM Pl )^^u(p 3 )(-tg s T a r)v(p 4 ) (5) 
5M r tex = ^(p 2 )(-^T a f v ) M ( Pl )^|^«(p 3 )(-^ s TV)^(P4) (6) 

^vertex = ^p 2 ) (-^ s T a 7 >(Pl) = ^u(Ps) H<? s T a f>(p 4 ) (7) 

s 

T 2 

+u(p 3 )P R u(p 1 )v(p 2 )PLv(p4)f2 + u{p 3 )P L u{pi)v{p 2 )PRv{p A )f- i 



+u{p 3 )P L u(p 1 )v(p 2 )P L v{p A )f 4: + u{p z )P R u{p 1 )v{p 2 )PR $3v(p 4 )f 5 

+u(p 3 )P R u(p 1 )v(p 2 )P L fov(p4)f 6 + u(p 3 )P L u(p 1 )v(p 2 )P R p , 3v(p4)f 7 

+u(p 3 )P L u(p 1 )v(p 2 )P L p , 3v(p4)f8 + u(P3) ^4Pru(Pi)v(p 2 )P R v(p 4 ) fg 
+u(Ps) p , APRu{pi)v(p 2 )P L v(p A )f 1 Q+ u(j) 3 ) f , APLu(pi)v(p 2 )P R v(p A )f 1 i 
+U(p 3 ) faPLU{pi)v{p 2 )PLV{pi)f 12 + U(p 3 ) jfAPRu{pi)v{p 2 )PR ^3^(P4)/l3 
+U(P3) P / 4Pru(Pi)v(P2)Pl P , 3v(p4)fu + U(p 3 ) ^Plu{Pi)v{P2)Pr p / 3 v(p4)fl 5 
+u(Ps) i>±PLu{$^)v{ip 2 )P L ^ 3 v{pi)fi & +u{p 3 )^^P R u{pi)v{p 2 )P R ^v{p A )fi 7 
+HP3h^PR u (Pl)^(P2)PLl^v(P4)fl8 + HP3h^ P Lu(pi)v(p 2 )P R ^v(p 4 )f 19 

+u{p 3 ) lll P L u{p l )v{p 2 )P L ^v{p A )f 2Q } (8) 



with 



= YP[ + Yl 5 P2 + k u F 3 + k^ 5 F' 4 + i k^F' 5 + i k^ l5 F^ (9) 
f" = 7 "Fi + -f-y 5 F 2 + k»F 3 + F75F4 + ik^Fs + ik^^Fe (10) 

Here, A; = pi + p 2 ,Pi,P2 denote the momenta of incoming partons,and p 3 ,p A are used 
for outgoing top-quark and its antiparticle. F-,Fi and fa are form factors which are 
presented in Appendix. 

The renormalized differential cross section of the subprocess is given by 

dcr 1 I3 t 2 



dcosO 327T s 



^|M ren | 2 , (11) 



where 



with 



|M ren | 2 = £|M | 2 + 2Re]TM + (5M] / + 5M 2 V + 5M b ), (12) 



£|M | 2 = g| [25m 2 + {i _ m 2 )2 + (fi _ m 2 )2]; (13) 
l> + <^ = ^§[2sm? + (f - m 2 ) 2 + (« - m 2 ) 2 ]^ - 1), (14) 
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5> + 5M 2 ^ = ||{[2Sm 2 + (t - m 2 ) 2 + (u - m 2 ) 2 ]^ - 1) + 2m 2 s 2 F 5 }, (15) 



£M+5M 6 = ^^{[sm 2 + (i-m 2 ) 2 ](f 2 + f 3 ) 

_ m ^2 + g _ 2m 2 s _ Q _ + £ _ m 2 )](/g + h _ fw _ /n) 

-[mis + (s - m 2 )(s + i- m 2 ) 2 + 2m 2 t (i - m 2 ) 

(s + t - m 2 )](f u + / 15 ) - 2{m 2 t s + (s + t- m 2 ) 2 ](f l8 + / 19 )}, (16) 

In the above equations, 9 is the scattering angle between the quark and the top quark, 
s, t, and u are the kinematic invariants for the 2—^2 subprocess with s + t + u = 2m 2 , 
and P t = \Jl — Amf/s the velocity of the final quarks. The total cross section for the 
production of top quark pair can be written in the form 

°~( s ) = £ / dx 1 dx 2 a ij (x 1 x 2 s,m 2 , ^i 2 )[F A (x l) fi)Ff(x 2 , n) + {A <->• 5)], (17) 

with 

s = (Pi + P 2 ) 2 , s = a;ia;2S, , . 

Pi = xxPx, p 2 = x 2 P 2 , (18) 

where A and B denote the incident hadrons, Pi and P 2 are their four-momenta, i,j 

are the initial partons, X\ and x 2 are their longitudinal momentum fractions, and the 

functions P/ 4 , Ff are the parton distributions of the initial-state hadrons A and B. In 

our numerical calculations, we have used the MRS Set A' parton distribution functions 

[8], and do not consider SUSY corrections to the parton distribution functions since 

the principle of decoupling demands that the corrections are negligible. Introducing a 

convenient variable r = xix 2 and changing to X\ and r as independent variables, the 

total cross section expression becomes 

^)=£/y(^)(^), (19) 

where r = 4m t /s and the quantity dL^/dr is the parton luminosity, which is defined 

as 



dLij f dxi 



<h J -^-F:H.r.. ll ).F! ! (j-., l ) + (A^ Pi. (20) 



T 



Now we present some numerical examples. In our numerical calculation, we input 
m t = 176 GeV and 1-loop a s (Q 2 = s) and use the phase space cuts: \rj\ < 2.5, pr > 
20GeV. As to the supersymmetric parameters involved in our calculations, in general, 
once tan (3 and mi L are fixed, we are free to choose two independent parameters in the 
stop mass matrix : (mf ,A t + //cot/3), which also can be tranfered to (mf , mfj. As 
explained above, we neglect the mixing of squark masses in the corrections to the gqq 
couplings and also neglect the mass splittings between squarks of different flavor for 
simplicity. 

In Figs. 1-3 we give some numerical results in a simple case, in which we set tan (3=1 
and mi L = m^ R = = rriq and thus the mixing angle is 7r/4 and A t + /i — m t . To 
compare the results in no-mixing case with those in the mixing case, we show the results 
in both cases in Figs. 1-3. From these figures we can see that the corrections in mixing 
case are smaller than in the no-mixing case. The plots versus squark and gluino masses 
in the mixing case get smoother than in no- mixing case. Fig.l show the dependence of 
the corrections on gluino mass for fixed squark mass of 150 GeV. The corrections can be 
either positive or negative, depending on the gluino mass. The corrections are negative 
for gluino mass below 150 GeV and above that they become positive. When gluino 
mass is changed from 100 GeV to 200 GeV, the corrections vary from -2%(-10%) to 
16%(23%) in mixing (no-mixing) case. The corrections get their positive maximum size 
at gluino mass of 200 GeV. When gluino mass is larger than 200 GeV, the corrections in 
both cases drop monotonously with the increase of gluino mass and tend to zero at 600 
GeV, showing the decoupling effects. The recent CDF lower limit [9] on gluino mass is 
160 GeV for arbitrary squark mass and 220 GeV when gluino mass is equal to squark 
mass. So the corrections are always positive if we consider the CDF limit on gluino 
mass. Fig.2 and Fig. 3 show the dependence of the corrections on squark mass for fixed 
gluino masses of 150 GeV and 200 GeV, respectively. The corrections in the mixing case 
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differ significantly from those in no-mixing case for low squark mass. For gluino mass of 
200 GeV the CDF lower limit [9] for squark mass is about 220 GeV. At this lower limit 
the corrections are 14% and 18% for the mixing and no-mixing cases, respectively. But 
with the increase of squark mass the difference of the corrections in both cases is getting 
negligibly small. 

From Figl-3 we have found that the corrections vary rapidly with gluino mass, es- 
pecially for 150 Gev < m~ g < 200 GeV, though mixied cases are smaller than unmixied 
ones. This is due to the fact that we have set m t = 176GeV in the numerical calculation, 
and the threshold for open top quark production is crossed in that region. If we change 
the top quark mass, we can find that such region is also shifted correspondingly, which 
provides a check on our calculation, especially on the treatment of the threshold. 

We also perform the numerical calculations for tan/? = 10. We found that the 
corrections are not sensitive to tan (3 value. For example, with gluino mass of 200 GeV 
and squark mass of 150 Gev we get Aa = 0.482 pb for tan/3 = 1 and Act = 0.490 pb for 
tan (3 = 10. 

In conclusion, we have shown that the supersymmetric QCD corrections to top quark 
pair production by qq annihilation in pp collisions can exceed 10% in both the mixing 
and no-mixing cases of squark masses even if we consider the recent CDF limit on squark 
and gluino masses. 

We thank J. Ellis for suggesting this calculations. This work was supported in part 
by the National Natural Science Foundation of China. 

Appendix 

The form factor Fi is given by 

Fl = l + %± \A U (F[ tfu) + 2m 2 t G ( ( fu) ) 

37T L V / 



+2m t m~ g B ii G { f i) + ^(F?> + 9Fj 



(i) 



<2h 



where 



An — CLiCLi + bibi, Bn — CLiCLi b{bi, 

ai = — b 2 = — ^=(cos# — sin^), a 2 = &i = — ^=(cos# + sin#), 
v 2 v 2 



JO 

= - I' dy 
Jo 



-m\y{\ -y) + m\{\ - y) + m|y n 



y n+1 (i~y) 



-mfy(l -y) +m|(l -y) + m 2 i y ) 



(21) 

(22) 
(23) 

(24) 

(25) 



= 2mi[m 5J B ii (c + c n ) - mtA^cu + c 2 i)] - 2AjjC 24 (-p 3 , k, m g , m it ,m k ), (26) 



F[ 2) = -1 + Aii[2c 2 4 + s(c 22 - c 23 ) - m?c - m 2 (c + 2cn + c 2 i)](-p 3 , mi v m~ g , m g 



-2B u m t mg(co + cn)(-p 3 , k,m,f.,mg,m,g), 



(27) 



Here 6* and 9' are the mixing angels of squarks. c , Qj are the 3-point Feynman integrals, 
definition and expression for which can be found in Ref. [10]. The form factor F[ is 
obtained by 



a. 



K j:j Fi qMj) + -(F± iy + 9F± 
8 



(i)' 



,(2)S 



where 



j _ A. jj \o-,e', Ljj - B JJ} 



F| = -2i^jQ24(pi,-A;,7n 5 ,77i^.,m^.), 



F x (2) = -1 + Kjj[2c24: + s(c 22 - c 23 ) - m~c ](pi, -A;, m^, m 5 , m 5 ) 



F 5 is given by 



24tt 



where 



[m t Aj(cn + c 2 i) - rrigBi^co + Cn)](-p 3 , fc, m 5 , m t ~ , m t ~. 



F 5 (2) = [m t Aj(cn + c 2i ) + m~gB ii c 11 \{-p 2l ,k,rni v m~g,m g: ), 



(28) 

(29) 
(30) 
(31) 

(32) 

(33) 
(34) 



fi are given by 

f 2 = m\o%D Q - m t m~ g o i3 o' l3 {D 12 + D 13 - D - D u ) + m\a%{D 2% - D 12 - D 24 ), (35) 

/ 6 = -mgCTijCr'ijiDo + D u ) + m t a^(D 12 + £> 24 ), (36) 
/io = m-gffij^D^ - D l2 ) - m t o%[D X2 + L> 23 + D 2 4 - ^13 - #26 - #25), (37) 

fu = al(D 12 + D 2A -D 13 -D 25 ), (38) 
As = 4^27, (39) 
(/3, h, fll, fl5, /19) = (/2, /6, /lO, /l4, /l8) aij-^Xij (40) 

where 

o-y = (a; + &i)(cj + dj), 4 = ( a * ~~ & *)( c i + d i)> 
Aij = (aj - bi)(cj - dj), X' i:j = (en + bi)(cj - dj), 

(c j ,d j ) = (a jf b j )\e^0 l ) (41) 
Here D , D i j(—p 1 , —p 2 ,p i ,m g ,mg j ,m g ,m^) are the 4-point Feynman integrals [10]. 
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Figure Captions 

Fig. 1 Relative correction to hadronic cross section at Tevatron versus gluino mass. 
Fig. 2 Same as Fig.l, but versus squark mass for gluino mass of 150 GeV. 
Fig. 3 Same as Fig. 2, but for gluino mass of 200 GeV. 
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